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ABSTRACT 

We apply our recently elaborated, powerful numerical approach to the high-resolution modeling of 
the structure and emission of circumstellar dust disks, incorporating all relevant physical processes. 
Specifically, we examine the resonant structure of a dusty disk induced by the presence of one planet. 
It is shown that the planet, via resonances and gravitational scattering, produces (1) an asymmetric 
resonant dust belt with one or more clumps intermittent with one or a few off-center cavities; and 
(2) a central cavity void of dust. These features can serve as indicators of a planet embedded in the 
circumstellar dust disk and, moreover, can be used to determine its major orbital parameters and even 
the mass of the planet. The results of our study reveal a remarkable similarity with various types of 
highly asymmetric circumstellar disks observed with the James Clerk Maxwell Telescope around Epsilon 
Eridani and Vega. The proposed interpretation of the clumps in those disks as being resonant patterns 
is testable - it predicts the asymmetric design around the star to revolve, viz., by 1.2°-1.6° yr~^ about 
Vega and 0.6°-0.8° yr"! about e Eri. 



Subject headings: circmstellar matter - dust, extinction 
: stars: individual (Vega, e Eridani 

1. INTRODUCTION 

A significant limitation to unambiguous planet detec- 
tion in circumstellar disks is the contaminating thermal 
emission from dust in the target systems. We propose to 
turn this hazard into an advantage for detecting and char- 
acterizing planetary systems. The major goal of the present 
Letter is to explore how the presence of the planet (s) im- 
pacts the disk via particular resonances responsible for the 
specific asymmetric features in the dust, and how these 
features would enable us to derive the major parameters 
of the yet invisible planetary system. The major obsta- 
cle in providing reliable, high-resolution numerical simu- 
lations - the particle-number limitation - is solved by our 
recently elaborated very efficient approach to numerical 
modeling of distributions of test particles in an external 
gravitational field (Ozernoy, Gorkavyi, & Taidakova 1998, 
2000; Gorkavyi et al. 2000; Taidakova & Gorkavyi 1999). 

2. RESONANCE STRUCTURES IN DUSTY 
CIRCUMSTELLAR DISKS 

Both the Poynting-Robertson (P-R) and stellar wind 
drags tend to induce dust inflow toward the star. As the 
dust passes by the planets in its infall, it interacts with 
them by accumulating in the outer planetary resonances, 
which, as we demonstrate below, could serve as an effi- 
cient means of planet detection. The existence of a reso- 
nance ring associated with the Earth's orbit was predicted 
by Jackson & Zook (1989); such a ring was indeed dis- 
covered in the IRAS (Dermott et al. 1994) and CORE 
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(Reach et al. 1995) data. The resonant design in the 
dust depends upon the capture probability, the particle 
lifetime in the resonance, the amplitude of libration, and 
the typical eccentricity of resonant particles. The total 
mass of dust contained in a particular resonance is the 
product of the probability for the dust particle being cap- 
tured into the resonance, Wf, the lifetime spent in the 
resonance, Tj (in a sufficiently dense environment, this 
time is due to particle collisions); and the total dust flow 
per unit time in the vicinity of the resonance, dM/At, i.e. 
Mj w WjTj{(lM/(it). Having introduced, for convenience, 
the mass of 'background' dust, « (1 - Wj)T(dM/dt), 
one gets M^/Mb w Wj{l-Wjy^ (tj/t), where r is the 
minimal timescale of the P-R drag (rpu), gravitational 
scattering (rgs), and collisional (tcou) characteristic times. 
Three particular cases are well distinguishable: 

1: T w Tpfl < Tgs,Tcoii- Actually, this coUisionless case 
is realized, for example, for Neptune in the Solar system. 
Assuming that Wj ^ 1 — Wj, one gets Mj/Mt ^ Tj/rpR. 
To get a discernable resonant structure, one needs to have 

Tj ^ TpR. 

2: r « Tgs < TcoU,TpR. This case is realized for a 
rather massive planet like Jupiter. If collisions are essen- 
tial (as they are for Vega's disk, see Backman 1998), then 
Tj ~ Tcou and one gets Mj/Mb ~ Tcou/TgaZl- The reso- 
nant structure has a high-contrast if Tgs <C TcoU- 

3: T « Tcoii < Tgs,Tpp. In this collisional case one has 
Tj « TcoU giving Mj/Mb ~ Wj/{1 — V^j), which results in 
a high-contrast resonant structure whenever Wj ~ 1. 
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It is worth emphasizing that if coUisions of resonant par- 
ticles arc important (cases 2 and 3), the resonant pattern 
can nevertheless exist (cf. Wyatt et al. 1999) as long as a 
source of particle replenishment is available. 

3. MODELLING CIRCUMSTELLAR DUST DISKS 
3.1. The Procedure 

The simplest planetary system under consideration, con- 
sisting of just one planet (orbiting a star of mass M at 
a circular orbit of radiiis a) possesses only mean motion 
resonances, withoiit most of secular resonances (except 
the Kozai resonances). We have computed about 300 
model disks in order to explore the effects of various res- 
onances adopting the following parameters: (1) the mass 
of the planet, m.pi: (2) the particle size coupled with stel- 
lar luminosity and mass, as described by the parameter 
/3 f« 0.3 -L*/(M*r), which is the ratio of the stellar light 
pressure and the gravitational force applied to the dust 
grain of density 2 g cm~^ and radius r (in /im), with L^, 
and Mi, the star's luminosity and mass in solar units; and 
(3) lifetime of dust particles, r. Below, we describe some 
of these models. 

3.2. Typical Resonant Structures Exemplified by the 
Resonances 1:1, 2:1, 3:1, and 3:2 

In the absence of a planet, the surface density of dust 
between the star and the pericenter of dust sources is con- 
stant, while the number density of dust varies as r~^, r 
being heliocentric distance (Gorkavyi et al. 1997). The 
presence of a planet induces dramatic changes in the dust 
distribution. This is illustrated in Fig. 1, which shows the 
geometry of dust distribution in different mean motion 
resonances with the planet. 

In the resonance 1:1, the orbits responsible for the dust 
design are of three kinds: (1) a horseshoe-like orbit that 
embraces the two Lagrange points L4 and L5 (see Fig. la); 
(2) an orbit with a small libration around L4 (Fig. lb, up- 
per feature): and (3) an orbit with a substantial libration 
around L5 (Fig. lb, lower feature). 

The dust structure in the resonance 2:1 looks similarly: 
the inner orbital loop precesses relative to L4 and L5 in the 
same way as the epicycle of particles in the 1:1 resonance. 
That loop can form a horse- like arc along the planet's or- 
bit (Fig. Ic) or can librate near either L4 or L5 (Fig. Id 
illustrates the case of a large libration). 

In the resonance 3:1, the dust particles librate around 
the L4 with a smaller amplitude (Fig. le); no horse- like 
libration has been observed. Finally, Fig. If illustrates 
the resonance 3:2 with an appreciable libration, which re- 
sults in a symmetric four-clump pattern. When libration 
is small, the resonance contains a two-clump pattern. 

In all cases mentioned above, the vicinity of the planet is 
free of dust. Each resonant dust structure is accompanied 
by other specific caverns in dust distribution (see Fig. 1). 

3.3. Dependence of Structure on Location of Dust Sources 

The scattered cometary population produces dust whose 
distribution in orbital elements tends to be rather smooth. 
A well-pronounced resonant structure in the dust can only 
be formed if there is a planet that massive enough to cap- 



ture some dust particles into resonances and to eject the re- 
maining non-resonant dust particles out (Case 2 in Sec. 2). 

Dust sources like the Kuiper belt objects (KBOs; here- 
after 'kuiperoids'), of which a substantial part is in res- 
onances, can produce the contrast resonant structures if 
(1) the drifting dust particles have a long lifetime and 
eventually are captured into resonances (the situation is 
similar to the formation of a resonant dust belt near the 
Neptunian orbit in the Solar system) or (2) both /3 and 
the lifetime of dust particles are small enough, so that the 
particles produced by the resonant kuiperoids cannot leave 
the resonances and thereby they form resonant structures 
(Case 3 in Sec. 2). Orbits of resonant objects look similar 
to dust orbits shown in Fig. 1. 

3.4. Dependence of Structure on Planetary Mass 

This dependence can be revealed from a variety of avail- 
able runs with different parameters. A low-mass planet, 
like Earth or Neptune, forms a resonant ring along the 
planetary orbit associated with high-order resonances, like 
3:2, 4:3, etc. (see, e.g. Fig. If). For more massive plan- 
ets, the resonances like 2:1, 3:1, etc. become essential, 
whereas the innermost resonances overlap (e.g. Wisdom 
1980). Since the hbration amplitude decreases as nipi in- 
creases, clumps similar to those seen in Fig.le (and not 
arcs like those shown in Fig. Id) would be typical for rUpi 
as large as mj, where mj is Jupiter's mass. 

4. A LINK TO OBSERVATIONS: IMAGERY OF 
CIRCUMSTELLAR DISKS 

Let us considcir some; observational signatures of planets 
embedded in disks. Figures 2a and 2b show the thermal 
emission from the simulated disks (seen face-on) with an 
embedded planet, and these are to be compared with the 
available observational data on Vega and e Eri shown in 
Figures 2c and 2d, respectively (the parent stars are seen 
almost pole-on). 

As a hotter star, Vega might result in a more exten- 
sive, compared to e Eri, dust production from evaporating 
comets. We adopt /3 = 0.3 (close to an upper limit of 
13 = 0.5 at which the dust particle escapes) and TcoU = 
2 X 10'' planet revolutions consistent with TcoU ~ 10^ yrs 
(Backman 1998). A planet as massive as ~ 2 mj induces 
the resonances n : 1, which results in two dusty clumps 
on either side of the planet, while the planet ejects the re- 
maining non-resonant dust particles out (Case 2 in Sec. 2). 
Thus, even if the distribution of the cometary population 
around Vega might be symmetric, ejection of non-resonant 
particles by the planet could make the resonant dust pat- 
tern highly asymmetric (Fig. 2a), which is reminiscent of 
the observed two asymmetric clumps near Vega shown in 
Fig. 2c (Holland et al. 1998). 

In contrast to Vega, e Eri is a less massive and a colder 
star with L « ^vega- Besides, the typical size of dust 
particles is several times larger than in Vega's disk (Back- 
man 1998), therefore we adopt (3 = 0.002 for the dust in 
e Eri's disk. Assuming that kuiperoids serve as the main 
source of dust in this disk, we find that the dust distribu- 
tion might reflect the distribution of the sources located 
mostly in the resonances such as 2:1 and 3:2. In this case, 
a comparatively small planetary mass of ^ 0.2 mj is able 
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Fig. 1. — Simulated surface density of the circumstellar dust (/3 = 0.065) captured into particular mean motion resonances with an outer 
planet shown as a square. The star's location is (0,0). In each panel, the surface density of the disk is shown on a linear contrast scale; the 
area covered is 200x200 (scale is arbitrary but, for the Solar system, it would be in AU); and the (single) planet is assumed to be on a circular 
orbit (ep;=0) at r = 30. The particle lifetime in the resonance, Tres, was taken to be 10'^ planet revolutions everywhere: a m = 0.3 mj, 
resonance 1:1 (a horscshoe-like orbit around the two Lagrangian points L4 and L5); 6 same as a, but an orbit around L4 with a minimal 
libration and another, nonoverlapping orbit around L5 with a maximum libratioii; c same as a, but resonance 2:1; d m = 0.3 mj, resonance 
2:1 (an orbit around L4 or L5 with a maximum libration); e m = 3 mj, resonance 3:1 (an orbit around the Lagrangian point L4 with a 
smaller libration amplitude caused by a larger mpi). Note that at rUpi as small as m = 0.3 mj, there is no capture into this resonance; / 
m = 0.3 mj, resonance 3:2. 



to induce, via libration, two arcs, with a pair of clumps 
at their edges (Fig. 2b). Even if the dust is distributed 
in equal proportions to both resonances, one could easily 
get an asymmetry seen in Fig. 2b by assuming that 10% 



of comets in the resonance 2:1 have a smaller libration 
amplitude in the lower branch of the pattern. A similar 
asymmetric structure has been revealed in the submillime- 
ter imagery of the e Eri disk shown in Fig. 2d (Greaves et 
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Fig. 2. — Thermal emission (A = 850 fim) from simulated circumstellar disks with the resonant structure dominated by a few resonances 
vs. observations. The planet shown in as a square in a, b revolves, like the disk, counterclockwise. In c, d, a square indicates the expected 
planet's location. The star's location is (0,0). To compare our simulations with observations, averaging with an appropriate beam size is 
applied to the simulations, a Simulation reminiscent of the disk of Vega: = 2.5 Mq; rUpi = 2 mj; (3 = 0.3; t = 2 X 10* planet's revolutions; 
the number of dust sources (comets) is 486, and the number of particle positions is 3 X 10^". b. The simulation reminiscent of the disk of 
e Eri: Mi, = 0.8 Mq; rripi = 0.2 mj; the resonances 2:1 and 3:2 (in equal proportions); (3 = 0.002; r = 10^ planet's revolutions; the number 
of dust sources (kuiperoids) is 3567, and the number of particle positions is 2 X 10^. c. Circumstellar disk of Vega (Holland et al. 1998). d. 
Circumstellar disk of e Eridani (Greaves ct al. 1998). 



al. 1998). 

Our modeling, which indicates that Vega may have a 
planet of mass ^ 2 mj at a distance of 50 — 60 AU, and 
e Eri may have a less massive planet of m w 0.2 mj at 
a similar distance of 55 — 65 AU, is testable. Each reso- 
nant feature is stationary in the reference frame corotating 
with the planet, but it is not so for the observer at Earth. 
Therefore, if our interpretation of asymmetric clumps as a 
dynamical resonant structure is correct, the above asym- 
metric feature revolves around the star (probably, coun- 
terclockwise) with an angular velocity of 1°.2-1°.6 yr~^ 
(Vega) and 0°.6-0°.8 yr~^ (e Eri) - a prediction that can 
be tested within several years. If confirmed, the proposed 



interpretation of the structure in Vega- and e Eri-like cir- 
cumstellar disks seen face-on would make it possible not 
just to reveal the embedded planet and determine its semi- 
major axis, but also to constrain its other basic parame- 
ters, such as the planet's mass and even to pinpoint the 
position of the planet (see Figs. 2a and 2b). 

5. CONCLUSIONS 

As we demonstrated above, highly asymmetric struc- 
tures in the circumstellar disk and a central 'hole' void 
of dust can serve as indicators of at least one planet em- 
bedded in the disk. The numerical simulations described 
above lead to the following major conclusions: 
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1 . The outermost planet in an exo-planetary system can 
produce: (a) an asymmetric resonant dust belt consisting 
of arc(s) and one or more clumps intermittent with one or 
a few cavities; and (&) a central cavity void of dust. 

2. The morphology of these belts and cavities (size, 
asymmetry, number of clumps, and their pattern) depends 
upon particular resonances involved and is eventually a 
function of the planetary mass, location of dust sources, 
and (3. 

3. The crucial test of the above picture would be the dis- 
covery of revolution of the resonant asymmetric structure 
around the star. For circumstellar disks in Vega and eEri 
the asymmetric design is expected to revolve, respectively, 
by 1°.2-1°.6 and 0°.6-0°.8 annually 

In sum, our numerical simulations offer a possibility to 
reveal the presence of the planet (s) in a circumstellar dusty 
disk seen pole-on and, under further development, could 
provide a method for determination of planetary parame- 
ters using the visible morphology of an outer part of the 
disk. The proposed approach seems to be sufficiently pow- 
erful to provide also the interpretation of disks seen edge- 
on in scattered light, such as the (3 Pic disk revealed by 
STIS optical observations (Heap et al. 2000). 
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